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1.  INTRODUCTION.  ' 

Let  us  consider  systems 

x  =  f(x,y,t,e)  , 

(1) 

cy  *  g(x,y,t,£)  , 

of  m  +  n  ordinary  differential  equations  on  a  finite  interval, 
say  0  <  t  <  1,  subject  to  q  initial  conditions 

A(x(0),y(0) ,e)  =  0  ,  (23 

and  r  terminal  conditions 

B (x(l) ,y (1) ,e>  =  0  ,  (3) 

with  q  +  r  =  m  +  n.  We  shall  assume  that  f,  g.  A,  and  B  have 
asymptotic  series  expansions  in  e  with  coefficients  being 
smooth  functions  of  the  remaining  variables,  and  we  shall  seek 
the  asymptotic  behavior  of  solutions  under  the  condition  that 
the  n  «  n  Jacobian  matrix  gy(x,y,t,0)  has  a  hyperbolic  split¬ 
ting  with  k  >  0  (strictly)  stable  and  n  -  k  >  0  (strictly)  un¬ 
stable  eigenvalues  for  all  x  and  y  and  for  0  <_  t  <_  1.  We 
shall  also  suppose  that  q  ».  k  and  r  ^  n  -  k,  since  linear 
examples  suggest  that  a  limiting  solution  as  the  small  posi¬ 
tive  parameter  e  tends  to  zero  is  unlikely  to  occur  otherwise. 

The  reader  should  realize  that  the  corresponding  initial 
value  problem  with  k  -  n  has  a  well-understood  asymptotic 
solution,  as  presented  in  Wasow  (1965)  and  O’Malley  (1974). 
Initial  value  problems  with  a  fixed  number  of  purely  imaginary 
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eigenvalues  have  been  considered  by  Hoppensteadt  and  Miranker 
(1976)  and  Kreiss  (1979),  while  Vasil 'eva  and  Butuzov  (1978) 
and  O'Malley  and  Flaherty  (1980)  discuss  problems  for  which  gy 
has  a  nullspaee.  We  note  that  such  problems  can  be  consider¬ 
ably  more  complicated  when  eigenvalues  of  gy  cross  or  approach 
the  imaginary  axis  (cf . ,  e.g.,  the  resonance  examples  of 
Ackerberg  and  O'Malley  (1970)  and  the  initial  value  problem  of 
Levinson  (1949)).  Finally,  note  that  the  strict  eigenvalue 
stability  assumptions  can  be  weakened  in  "boundary  layer 
regions"  (cf.  Howes  and  O'Malley  (1980)).  The  two-point  prob¬ 
lems  arise  naturally  in  optimal  control  theory  (cf.  Kokotovic 
et  al.  (1976)  and  O'Malley  (1978)),  among  many  other  applica¬ 
tions.  Moreover,  knowing  about  the  asymptotic  behavior  of 
solutions  is  extremely  helpful  in  developing  schemes  for  the 
numerical  solution  of  stiff  boundary  value  problems  (cf. 

Hemker  and  Miller  (1979)  and  Flaherty  and  O'Malley  (1980)). 

2 .  ' THE  ASYMPTOTIC  APPROXIMATIONS. 

With  the  assumed  hyperbolic  splitting,  we  must  expect  so¬ 
lutions  to  feature  nonuniform  convergence  as  t  -  0  (i.e., 
boundary  layers)  near  both  endpoints.  Indeed,  it  is  natural 
to  seek  bounded  (uniform)  asymptotic  solutions  in  the  form 


x(t,e)  "  X(t.c)  +  e5(t,e)  +  en  (o ,  e )  , 
y(t,c)  -  Y  ( t ,  c  )  ♦  u(t,e)  ♦  v(o,e)  , 


(4) 


°n  0  <  t  <  1,  where  the  outer  solution  (X(t,c) ,Y(t,E) )  repre¬ 
sents  the  solution  asymptotically  within  (0,1),  where  the 
initial  layer  correction  (e£  (r,c) ,u (t,e) )  decays  to  zero  expo¬ 
nentially  as  the  stretched  variable 

x  »  t/c  (5) 

tends  to  infinity,  and  where  the  terminal  layer  correction 
(cn (c , e ) , v (o , e ) )  goes  to  zero  as 

o  *  (1  -  t)/t  (6) 


becomes  infinite.  Within  (0,1),  then,  such  solutions  are  re¬ 
presented  by  an  outer  expansion 


X(t,C) 

Y(t.E) 


(7) 


I 


V 


.  m  jwyi 


Xj  (t) 
Yjlt) 


The  limiting  uniform  approximation  corresponding  to  (4)  is 


x(t ,£)  »  XQ(t)  +  0(e)  , 

y(t,£)  «  Yq (t)  +  Uq(t)  +  Vg(o)  +  0(e)  » 

on  0  <  t  <  1.  At  t  »  0,  the  singularly  perturbed  or  fast 
vector  y  usually  has  a  discontinuous  limit,  jumping  from 
y(0,0)  »  Yq ( 0)  +  Ug(0>  to  Yq(0)  at  t  -  0+.  An  analogous 
Heaviside  discontinuity  occurs  near  t  =  1  whenever  Ug(0)  0, 

and  the  derivative  y(t,£)  generally  features  delta-function 
type  impulses  as  e  *  0  at  both  endpoints.  (The  relation  of 
such  observations  to  linear  systems  theory  is  of  considerable 
current  interest  (cf.  Francis  (1979)  and  Verghese  (1978)).  For 
problems  linear  in  the  fast  variable  y,  we  can  also  find  un¬ 
bounded  solutions  with  endpoint  impulses  (cf.  Ferguson  (1975) 
and  the  Appendix  to  this  paper) . 

The  outer  expansion  (7)  must  satisfy  the  full  system  (1) 
within  (0,1)  as  a  power  series  in  e.  Thus,  the  limiting  solu¬ 
tion,  (Xq,Yq) ,  will  satisfy  the  nonlinear  and  nonstiff  reduced 
system 

XQ  *  f(Xo,Yg,t,0)  ,  0  »  g(X0,Yg,t,0)  (9) 

there.  Because  gy  is  nonsingular,  the  implicit  function 
theorem  guarantees  a  locally  unique  solution 


Y0(t)  -  G(X0,t) 

of  the  latter  algebraic  system,  so  there  remains  an  m-th  order 
nonlinear  system 

Xg  -  F(Xg,t)  S  f  (Xg,G(Xg,t)  ,t,0)  (ll) 

for  Xg.  Later  terms  of  the  expansion  (7)  will  satisfy 
linearized  versions  of  the  reduced  system.  The  coefficients 
of  c  provide  that 

*1  "  fxtX0'Y0't,0) h  +  fy(X0'Y0’t,0,Yl  +  £E(X0'Y0't,0)  ' 

and 

Y0  ’  gx(X0'Y0't'0>Xl  *  gy(xo'Y0't,0)Yl  +  9e (X0'Y0't,0>  ' 

so  Yl(t)  -  G^  ♦  g;1^  -  gE)  and  ^  -  F^  ♦  -  9C> 

+  f  .  More  generally,  for  each  X  >  1,  we’ll  obtain  a  system 
of  the  form 
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Ylc<t>  "  GJC(X0.t)X){(t>  t-  ok_1(X(),...,Xk_1.t)  , 
Xk(t>  “  Fx<X0't>Xk(t)  +  9k-l(X0*  ‘  " 'Xk-l,t:)  ' 


(12) 


with  successively  determined  nonhomogeneous  terms. 

In  order  to  completely  specify  the  outer  expansion  (7), 
we  must  provide  boundary  conditions  for  the  m  vectors  Xk(t), 
k  >  0.  Most  critically,  we  first  need  to  provide  m  boundary 
conditions  for  the  "slow"  vector  XQ(t)  »  x(t,0)  in  order  to 
determine  the  limiting  solution  (Xq,Yq)  within  (0,1).  It  may 
be  natural  to  attempt  to  determine  them  by  somehow  selecting 
some  subset  of  m  combinations  of  the  m  +•  n  boundary  conditions 
(2)  and  (3)  evaluated  at  e  ■  0  (cf.  Flaherty  and  O'Malley 
(1980)  where  this  is  done  for  certain  quasilinear  problems). 
For  scalar  linear  differential  equations  of  higher  order,  the 
first  such  cancellation  law  was  obtained  in  Wasow's  NYU  thesis 
(cf.  Wasow  (1941,  1944)).  For  linear  systems  (with  coupled 
boundary  conditions) ,  a  (necessarily)  more  complicated  can¬ 
cellation  law  is  contained  in  Harris'  postdoctoral  efforts 
(cf.  Harris  (1960,  1973)).  These  significant  early  works 
suggest  that  we  should  seek  a  cancellation  law  which  ignores 
an  appropriate  combination  of  k  initial  conditions  and  of 
n-k  terminal  conditions,  so  the  limiting  solution  is  deter¬ 
mined  by  a  nonlinear  m-th  order  reduced  boundary  value  problem 


X0  -  F(XQ,t)  ,  0  <  t  <  1 

»(XQ(0))  -  0  ,  t (XQ (1) )  -  0  , 


(13) 


involving  q  -  k  initial  conditions  and  r  -  n  ♦  k  terminal  con¬ 
ditions.  Hoppensteadt  (1971)  considered  the  reverse  problem: 
Given  some  solution  of  a  reduced  problem,  what  conditions 
guarantee  that  it  provides  a  limiting  solution  to  the  original 
problem  (l)-(3)  within  (0,1).  Hadlock  (1973)  and  Freedman  and 
Kaplan  (1976)  also  consider  singular  perturbations  of  a  given 
reduced  solution,  as  do  Sacker  and  Sell  (1979)  who  allow  the 
Jacobian  g^  to  be  singular  but  subject  to  a  "three-band" 
condition.  Vasil'eva  and  Butusov  (1973)  consider  problems 
with  special  boundary  conditions,  though  their  results  are  ex¬ 
tended  to  more  general  boundary  conditions  by  Esipova  (1975). 
We  note  that  the  numerical  solution  of  a  reduced  problem  like 


(13)  is  much  simpler  than  that  of  the  original  problem, 
because  (13)  is  not  stiff  and  its  order  is  n  instead  of  n  t  > 
Thus,  the  solution  of  (13)  might  be  used  as  an  approximate  so 
lution  of  the  original  (full)  problem  from  which  to  obtain 
better  approximations  by  adding  boundary  layers  and  by  using 
Newton's  method. 

Near  t  ■  0,  the  terminal  boundary  layer  correction  is 
negligible,  so  the  representation  (4)  of  our  asymptotic  solu¬ 
tions  requires  the  initial  layer  correction  (e£,u>  to  satisfy 
the  nonlinear  system 


d£  dx  _  dx 
li  3t  3t 


f ( X  +  e£,Y  *  u.ct.c)  -  f(X,Y,ex,e)  , 


a?  *  e(at  -  as>  ’  9(x  +  ee'Y  +  u'eT'G)  -  9(x'Y'eT'e)  ' 

on  r  >  0  and  to  decay  to  zero  as  i  -  *>.  This,  in  turn,  pro¬ 
vides  successive  differential  equations  for  the  coefficients 
in  the  asymptotic  expansion 

f£(T,E)]  -  i 

-I  3  .  (15 

u(r, e)J  j —0 j  ( t ) 

Thus,  when  e  »  0,  we  have  the  limiting  initial  layer  system 


-  f (X0(0),Y010)  +  yQ, 0,0)  -  f (XQ (0) , Y„ (0) ,0,0)  , 


-g~-  -  g (XQ ( 0) , Yg ( 0)  +  u0,0,0)  -  g(X„(0),Y0(0) ,0.0) 
The  decay  requirement  determines 


as  a  functional  of  uQ,  while  uQ  satisfies  a  conditionally 
stable  nonlinear  system 

duo 

-  g(Xjj(0)  ,Yq(0>  +  U0,0,0)  5  G0(u0:x(0,0)  )u0  .  (17 


Our  hyperbolicity  assumption  on  the  eigenvalues  of  gy  there¬ 
fore  implies  that  the  limiting  boundary  layer  correction  is 
determined  by  a  classical  conditional  stability  problem  (17) 


on  t  >  0.  The  standard  theory  (cf.  Coddington  and  Levinson 
(1955)  or  Hartman  (1964)  or,  in  more  geometrical  terms, 
Fenichel  (1979)  and  Hirsch  et^  a_l.  (1977))  shows  that  for  each 
x(0,0),  there  is  (at  least  locally)  a  k-manifold  I(x(0,0>) 
nontrivially  intersecting  a  neighborhood  of  the  origin  such 
that  for 

u0(0)  e  I (x ( 0, 0) )  (18) 

the  initial  value  problem  for  (17)  has  a  unique  solution  on 
t  >  0  which  decays  to  zero  exponentially  as  r  ♦  One  very 
difficult  problem  is  how  to  compute  the  stable  initial  mani¬ 
fold  I,  even  when  x(0,0)  is  known.  Hassard  (1979)  has  begun 
to  address  this  problem  through  a  Taylor's  series  approach  and 
Kelley's  representation  of  such  stable  manifolds  through  the 
center  manifold  theorem. 

Recalling  that  the  q  limiting  initial  conditions  take  the 

form . 

A ( XQ ( 0 ) , G ( Xg ( 0 ) ,0)  +  uQ(0),0)  ■  0  (19) 

(cf.  (2)  and  (8)),  we  will  assume  that  it  is  possible  to  solve 
k  of  these  q  equations  (perhaps  nonuniquely)  for  an  isolated 
solution 

uQ(0)  5  y(X0(0))6  I ( x0 (0) )  .  (20) 

Phrased  somewhat  differently,  in  the  style  of  Vasil 'eva  (1963), 
we  are  asking  that  the  initial  vector  u0<0)  for  the  leading 
term  of  this  initial  layer  correction  belong  to  the  "domain  of 
influence"  of  the  equilibrium  point  Uq(t)  “  0  °f  the  initial 
layer  system  (17)  which  is  itself  parameterized  by  x(0,0)  » 
Xg(0).  Rewriting  the  remaining  q  -  k  initial  conditions  as 

*(X0  <0) )  5  VA(X0 ( 0) , G (XQ ( 0) , 0)  *  y(X0(0U ,0)1 '  »  0  (21) 

(where  the  prime  indicates  the  appropriate  q  -  k  dimensional 
subvector) ,  we  thereby  specify  the  initial  conditions  needed 
for  the  reduced  boundary  value  problem  (13).  Because  (21) 
generally  depends  on  y,  we  note  that  the  conditions  (21)  do 
not  simply  correspond  to  a  subset  of  the  original  initial  con¬ 
ditions  evaluated  along  the  limiting  solution  (Xg,Yg). 

In  the  important  quasilinear  case  when  gy(x,y,t,0)  » 
g(x,t)  is  independent  of  y  (at  least  near  t  »  0) ,  the 


resulting  initial  layer  system  (17)  is  linear  and  wq(t)  - 

V 

e  u0<0)  for  Gg  *  Gq(Xq(0),0).  Thus,  uQ  will  decay  to  zero 
as  t  -  •  if  U0(0)  *  Pg (Xg ( 0) ) u g (0)  where  PQ  ■  Pg  projects  onto 
the  k  dimensional  stable  eigenspace  of  Gg.  (More  generally, 
the  manifold  I  will  not  coincide  with  the  stable  eigenspace  of 
Gg.)  when  we  further  assume  that  Ay ( XQ ( 0 > ,y,0)  I  A(Xg(0)>  is 
independent  of  y,  there  will  be  a  unique  solution  Ug(0)  in  the 
then  fixed  manifold  I(Xg(0)),  provided  the  matrix  A(Xg(01) 

• Pg (Xg ( 0))  has  full  rank  k  let.  Flaherty  and  O'Malley  (1980)). 

A  simple  nonlinear  example  occurs  when  y  is  a  scalar  and 
g  (x,y ,  t,  0)  =*  g1(x,t)y2  +  g2(x,t)y  +  g3(x,t).  Then  uQ  will 
satisfy  a  Riccati  equation  with  solution 

VT)  “  Vo(0,/l(Ho  +  Volo))e  H°T  "  Vo(on 

for  G-g  »  g^fXgfOJ.O)  and  »  g^  (XQ  (0)  ,  tQ  (0)  ,  0, 0)  as  long  as 
the  denominator  is  nonzero.  If  HQ  >_  0,  only  the  trivial 
initial  layer  correction  Ug(r)  =  0  will  be  zero  at  infinity. 
With  the  stability  assumption  HQ  <  0,  however,  existence  on 
r  0  and  exponential  decay  at  infinity  is  guaranteed  provided 
Hg  +  GgPg(O)  <  0-  Thus,  the  magnitude  of  the  initial  layer 
jump  Ug(0)  must  be  restricted  when  GgWg(0)  >  0. 

The  terminal  layer  correction  (en  (o,e) ,v(o,c) )  can  be 
analyzed  quite  analogously  to  the  initial  layer  correction.  In 
particular,  the  leading  terms  (hg.Vg)  will  be  determined 
through  exponentially  decaying  solutions  of  the  conditionally 
stable  terminal  layer  system 

dvo 

-3^  -  -g(Xg(l),Yg(l)  +  Vg,l,0>  =  -G^Vg.-xd.ODVg  (22) 

on  o  >  0,  which  has  an  n  -  k  dimensional  manifold  T(x(l,0))  of 
initial  values  Vg(0)  providing  decaying  solutions  to  (22)  as 
t  ■»  ».  If  we  then  assume  that  n  -  k  of  the  r  limiting 
terminal  conditions 

B(Xg(l),G(Xg(l) ,1)  +  vQ(0),0)  *  0  (23) 

provide  an  isolated  solution 

v0(0)  5  6 (Xg(l»  €  T(X0<1>)  ,  (24) 

the  remaining  r  -  n  +  k  conditions  provide  the  terminal 
conditions 


(25) 


! 

;;  :>i 
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fUgll))  =  (B(X0(1)  ,G(X0(1)  ,1)  +  «  (XQ  (1)  )  ,  0)  - 


for  a  reduced  boundary  value  problem  (13) . 

The  reduced  two-point  boundary  value  problem  (13)  con¬ 
sists  of  the  nonlinear  reduced  equation  (11)  of  order  m 
together  with  the  m  separated  nonlinear  boundary  conditions 
(21)  and  (25).  If  it  is  solvable,  such  a  reduced  problem  can 
have  many  solutions.  Corresponding  to  any  of  its  isolated 


solutions  Xq ( t) ,  one  can  expect  to  obtain  a  solution  of  the 


original  problem  (l)-(3),  for  e  sufficiently  small,  which  con¬ 
verges  to  (Xg,G(X0,t))  within  (0,1)  as  e  -  0.  Sufficient 
hypotheses  on  the  corresponding  linearized  problem  to  obtain  a 
uniform  asymptotic  expansion  (4)  are  provided  by  Hoppensteadt 
(1971)  and  others.  For  this  reason,  we  shall  merely  indicate 
the  considerations  involved  in  obtaining  further  terms  of  the 
initial  boundary  layer  correction  and  boundary  conditions  for 
higher  order  terms  of  the  outer  expansion. 


Further  terms  of  the  initial  layer  correction  (15)  are 

k  in  the 


determined  from  the  corresponding  coefficients  of  e 
nonlinear  system  (14).  Thus,  we  must  have 

dC 


sr  •  yvo)'V°> 


P  g  (  X  )  fO.OJU^  +  Pfc-ltf) 


duk 

■ar 


(26) 


yx0(°)'V°> 


+  U0  (T  )  ,0,0)uk  +  CJlc_i  (  T  ) 


for  k  >_  1,  where  the  nonhomogeneous  terms  will  be  exponen¬ 
tially  decaying  as  x  ♦  •  because  the  preceding  C^’s,  p^'s,  and 
their  derivatives  so  behave.  The  homogeneous  systems  are 
linearizations  of  that  for  (Cq,Uq),  end  the  decaying  vector  £k 
will  be  uniquely  provided  in  terms  of  uk  by 

“  d£. 


£k(x) 


-I 


dx 


(s)ds 


(27) 


To  obtain  p^,  it  is  natural  to  first  consider  the  variable 


coefficient  homogeneous  system 
g(x)u 


du 

dx 


(28) 


with  g (x )  -  gy ( XQ ( 0) ,VQ(0)  +  jg  (x ) ,  0,  0)  .  Our  hyperbolicity 
assumption  (more  specifically,  the  eigenvalue  split  for  g(») 


and  the  exponential  convergence  of  g(T)  to  g(»))  guarantees 
that  (28)  will  have  k  linearly  independent  exponentially 
decaying  solutions  js  t  *  •.  We  assume  that  the  split  is 
maintained  for  all  t  >  0.  If  we  let  U(t)  be  a  fundamental 
matrix  for  (28)  with  U(0>  =  I  and  let  PQ  be  a  constant  n  *  n 
matrix  of  rank  k  such  that  U(t)Pq  provides  the  linear  subspace 
of  decaying  solutions  to  (28),  the  decaying  solution  uk  of 
(26)  must  be  of  the  form 

uk(T)  =  U(T)PQck  +  uk-i(T>  (29) 

with  the  particular  solution 

T 

Wk-1(T)  *  {  U<T)P0o"1(r)qk_1(r)dr 

0 

CO 

-  |  U(t)(I  -  P0)u-1(r)qk_1(r)dr  . 

T 

The  vector  ck  remains  to  be  determined.  In  problems  where  GQ 
is  constant,  Pq  is  the  P^  used  for  the  quasilinear  problem. 

The  use  of  such  exponential  dichotomies  (cf .  Coppel  (1978) )  in 
the  singular  perturbations  context  goes  back  to  Levin  and 
Levinson  (1954).  Indeed,  the  "roughness"  of  the  exponential 
dichotomy  might  be  used  to  justify  the  use  of  U(t)pq  all  the 
way  back  to  r  *  0.  Proceeding  analogously,  the  terminal  layer 
term  vk(c)  will  be  determined  in  the  form 

vk(o)  =  VfcJPjdj^  +  v£_i_(°>  *  (30) 

where  V(a)P^  is  assumed  to  span  an  n  -  k  dimensional  space  of 
decaying  solutions  to 

a?  ’  -VV^'V1)  +  V°)'i'0>v 

on  a  >  0  and  uk_^(°)  i*  exponentially  decaying  and  success¬ 
ively  determined.  nk(o)  will  uniquely  follow  from  vk  by 
integrating  dnk/da.  To  complete  the  formal  determination  of 
our  expansion  (4),  we  must  successively  specify  the  constants 
c.  in  (29),  d.  in  (30),  and  the  m  boundary  conditions  for  each 
Xk(t) ,  k  >  0. 


Since  x(0,e)  -  X(0,e)  +  cC(0,c)  and  y(0,c)  '  Y(0,e) 

+  u(0,c),  the  coefficient  of  (for  any  k  >  0)  in  the  initial 
condition  (2)  implies  that 


W0>'V0)  +  lJo(0>'0)Xk<0) 

+  Ay(XQ(0)  ,V0(0)  +  u0  (0)  ,  0)  (Vk  (0)  +  uk(0)) 

is  successively  determined.  Since  V^IO)  -  Gx (XQ (0) , 0) ( 0) 
and  Mk<0)  -  are  also  known  in  terms  of  preceding 

coefficients,  we  have 


(Ax0  +  Ay0Gx0)Xk(0)  +  AyOP0ck 


determined  termwise.  (The  2ero  subscripts  indicate  evaluation 
along  <x(0, 0) ,y (0, 0) , 0) . )  Assuming  that  the  matrix 


A  <x(0,0),y(0,0),0>P0  (32) 

has  its  maximal  rank  k,  it  will  be  possible  to  (perhaps  non- 
uniquely)  solve  k  of  the  q  equations  (31)  for 

ck  «  P0=k  '  (Ay0P0)+t6k-l  -  (Ax0  +  Ay0Gx0)Xk(0>>  (33) 

(with  the  dagger  representing  the  matrix  pseudoinverse) .  This 
leaves  the  remaining  q  -  k  initial  conditions 

rxk(0>  =  Ck_x  (34) 

to  be  solved  for  Xk { 0 ) .  Here 

Y  -  TI  -  AyoP0(AyOP0)  1 <Ax0  +  Ay0Gx0>  ' 

ek-1  =  Ay0p0<AyoPo,+ik-l’  and  y  will  have  rank  q  -  k. 

In  analogous  fashion,  if  the  matrix 

By(x(l,0) ,y(l,0) ,0)P1  (35) 

has  its  maximal  rank  n  -  k,  we  can  use  n  -  k  of  the  terminal 
conditions  (3)  to  (generally  nonuniquely)  provide  d^  »  P.d^  i" 
(30),  and  the  remaining  r  -  n  +  k  terminal  conditions 

«Xk<l)  »  nk_1  (36) 


for  the  outer  expansion  term  X^lt).  Here,  nk_1 
successively  and 


is  known 


[I  -  BylPl(BylPl)](Bxl  ♦  BylGxl> 


with  the  subscript  1  indicating  evaluation  at 
(x  (1, 0) , y ( 1 , 0) ,1) . 

Putting  everything  together,  we've  shown  that  the  k-th 
term  in  the  outer  expansion  should  satisfy  an  m-th  order 
linear  boundary  value  problem 


Xk  =  FxfX0't)Xk  +  Sk-l(X0 . Xk-l'fc)  ' 

*V0)  -  5k_i  .  *VD  ■  Vl  ' 


(37) 


with  successively  determined  nonhomogeneities  6^  and 

These  problems  for  all  k  >  0  will  have  unique  solutions 
Xj^  ( t )  provided  the  corresponding  homogeneous  system 

X  *  F  ( X  » t)  X  , 

yX  { 0)  *  0  ,  6  X  ( 1)  =  0  , 


has  only  the  trivial  solution  on  0  <_  t  <  1.  Otherwise,  one 
must  impose  appropriate  orthogonality  conditions  on  (37)  to 
achieve  solvability. 

Altogether,  we've  determined  possibilities  for  con¬ 
structing  multiple  formal  asymptotic  solutions  to  our  boundary 
value  problem  ( 1 )  —  ( 3 )  in  the  form  (4).  To  prove  that  the 
corresponding  solutions  exist  for  e  sufficiently  small  re¬ 
quires  some  further  analysis  (cf.  Hoppensteadt  (1971), 
Vasil'eva  and  Butuzov  (1973),  and  Eckhaus  (1979)),  but  no 
surprises. 

3.  A  SIMPLE  EXAMPLE. 

A  relatively  simple  example  is  provided  by  the  harmonic 
oscillator  system 

x  =  1  ~  x  , 

eyx  =  y2  ,  (38) 

ty2  =  a2 (xly^  +  S (x)  , 

where 


a(x)  =  1  +  2x 


and 


8 (x)  =  8x(l  -  x)  . 


TV 
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Since  m  =  k  =  n-  k  =  l  whenever  a(x)  remains  nonzero,  we  can 
expect  one  dimensional  boundary  layer  behavior  at  each  end¬ 
point.  Using  the  three  linear  boundary  conditions 

x ( 0)  +  y^fO)  »  0  , 

-bx(0)  +  y2 ( 0)  =  0  ,  (39) 

x(l>  +  y1(l)  -  0  , 

we  could  hope  that  the  limiting  interior  behavior  is  deter¬ 
mined  by  the  reduced  system 


y20  -  0  ,  (40) 

u2(X0)Yio  +  S(XQ)  =  0 

and  a  combination  of  the  limiting  initial  forms  Xq ( 0 )  +  Y^Q(0) 
and  -bXQ(0)  set  to  zero. 

By  expJiciting  analyzing  the  initial  layer  system  and 
determining  the  appropriate  projection  matrix  P.,  we  find  the 
initial  condition  ♦ ( XQ ( 0 } )  *  0  (needed  for  the  reduced  solu¬ 
tion)  to  be 

i ot (X0 ( 0)  )  |  (X0 ( 0)  +  V1Q(0))  -  bXQ(0)  =*  0  .  (41) 

Rewriting  this  as  a  cubic  polynomial,  we  obtain  the  three 
roots 

X0 ( 0)  =  0 

and 

Xg(0)  *  -  j  ^  sgn  o0|b  +  Vtb  sgn  aQ  -  4) 2  +  48^ 

arranged  so  that  sgn  *  +1  depending  on  the  sign  of 
1  +  2Xg(0).  All  roots  are  appropriate  for  determining  an  XQ 
except  for  that  range  of  initial  values  for  which  the  re¬ 
sulting  a(Xg(t))  »  3  +  2(Xq(0)  -  l)e  t  has  a  zero  in  0  <  t  <  1. 
Then,  the  equation  for  y2  has  a  turning  point  and  Yjg(t) 
becomes  unbounded.  When  bXg(0)  =*  0,  both  of  the  initial 
conditions  for  the  system  are  satisfied  by  the  limiting 
solution,  and  there  is  no  nonuniform  convergence  at  t  *  0, 


i 


Otherwise,  there  is  an  initial  layer  with  a  nontrivial  Uq(t). 
Integrating  XQ ( t )  to  t  ■  1  (or  evaluating  the  explicit 
solution)  also  determines  Y^g(l).  If  Xg(l)  +  T^q(1)  ^  ®*  the 
terminal  layer  correction  Vg(o)  is  also  nontrivial. 

For  b  «  2,  we  obtain  the  three  roots  Xg(0)  «  0,  0.803, 
and  -4.29  and  three  corresponding  asymptotic  solutions.  For 
small  values  of  s,  the  asymptotic  solution  can  be  used  as  a 
first  guess  to  obtain  numerical  solutions  (cf.  Flaherty  and 
O'Malley  (1980)  and  Figures  1-3).  For  b  »  0  and  XQ ( 0 )  *  -7/2, 
a(XQ(t))  has  a  zero  above  t  *  1.  This  forces  the  boundary 
layer  jump  |y^q(0)  -  y  ( 1 , 0 )  |  to  be  large  (89.8  compared  to  0.4 
and  0.6  for  the  other  b  ■  0  roots,  Xg(0)  •  0  and  1/2)  (cf. 
Figure  4) . 

APPENDIX:  THE  CONSTRUCTION  OF  IMPULSIVE 
SOLUTIONS  TO  QUASlLINEAR  PROBLEMS. 

It  is  basic  to  the  preceding  development  that  we  can 
solve  the  limiting  boundary  conditions  for  initial  values  of 
the  boundary  layer  correction  terms  (i.e.  u Q ( 0 )  in  (19)  and 
vQ(0)  in  (23)).  Moreover,  the  solutions  must  lie  on  the 
stable  manifolds  I  and  T  for  the  corresponding  boundary  layer 
systems.  This  will  certainly  be  possible  when  the  limiting 
boundary  conditions  A (x (0) , y (0) , 0)  or  B (x (1) , y ( 1) , 0)  are  inde¬ 
pendent  of  y  and  suggests  that  the  corresponding  asymptotic 
solutions  will  then  no  longer  have  the  form  (4).  Indeed,  the 
earlier  work  of  O'Malley  (1970)  (with  m  «  n  »  1)  suggests  that 
the  endpoint  behavior  will  then  be  more  singular  (i.e., 
impulsive).  With  the  resulting  unboundedness  in  the  y  vector, 
we  ask  for  linearity  in  the  fast  variables  in  order  to 
generate  expansions. 

As  a  sample  problem,  consider 

3?  “  V*'6'*)  +  £2(t'c>*  ' 

dv  (42) 

e  St  *  +  g2(t,t)y  , 

where  g2<t,0)  has  k  >  0  stable  and  n  -  k  >  0  unstable  eigen¬ 
values  throughout  0  <  t  <  1,  subject  to  the  q  +  r  ■  m  +  n 
boundary  conditions 


'v*  . 


•*S*V«C 


I 


r 


A(x(0)  ,y  (0)  ,0  5  A^xfOl.c)  ♦  eA2(e)y(0)  -  0 
B ( x ( 1 ) ,y (1)  ,e)  »  0  . 


(43) 


Since  Ay  =  0  at  e 
tions  in  the  form 

x(t,c)  -  X( t, c )  ♦ 

y(t ,£)  ■  Y(t,c)  ♦ 


•  0,  we  shall  now  seek  asymptotic  solu- 


C (x,c)  +  en(o,e)  , 
i  y  (t  ,  e)  +  v(o,e)  , 


(44) 


(as  an  alternative  to  (4))  where  the  terms  are  all  expandable 
in  c  and  the  boundary  layer  corrections  decay  to  zero  as 
before . 

The  limiting  solution  (Xg,Yg)  within  (0,1)  satisfies  the 
reduced  system 

dx 

-  F(XQ,  t)  =  f^Xg.t.O)  4-  f  2  ( t ,  0 )  G  ( Xg ,  t)  ,  (4  5) 

where 

YQ(t)  -  G(X0,t)  S  -g^lt.OJg^Xg.t.O)  .  (46) 

The  initial  layer  correction  (C,i  u)  now  satisfies  the  almost 
linear  system 

^  »  f2(cx,c)u  +  £  [  f  ^  ( X  ♦  C,£T',E)  -  f^(X,CT,C)l  , 


“  g2(et,£)u  +  c(g1(X  +  C,tr,c)  -  gj^X.et.e)]  , 


so  its  leading  coefficients  satisfy  the  constant  linear  system 


du. 


F0y0  ' 


■ffT  ■  Gowo  ' 


for  (Fg.Gg)  »  ( f 2 (0, 0) ,g2 (0, 0) ) .  The  decaying  solutions  are 


U0(T>  "  *  °  Vo(0)  ' 
VT)  *  roGVo(T)  ' 


(47) 


where  PQ  projects  onto  the  constant  k  dimensional  stable 
eigenspace  I  of  Gg.  These  vectors  and  higher  order  terms  in 
the  initial  layer  correction  will  necessarily  lie  in  this  same 
(known)  eigenspace.  With  the  expansion  (44),  our  initial 
condition  takes  the  limiting  form 


I 


Ax(X0<0)  ♦  P0G0aP0u0(0>,0)  «•  A2(0)P0u0(0)  -  0  .  (48) 

w«  with  to  solve  k  of  these  q  nonlinear  equations  for 

VO101  5  Y(X0<0>)  • 

The  solution  will  be  locally  unique  if  the  appropriate  k  *  k 
Jacobian  is  nonzero.  The  remaining  q  -  k  initial  conditions 
will  then  provide  the  initial  values 

♦  <X0<0>>  5  (A(X0(0)  +  FgG'SUgtOn.O) 

(49) 

♦  a2 (0)Y (X0 (0) ) ) •  -  0  , 

for  the  limiting  solution  XQ(t). 

Preceding  analogously,  the  limiting  terminal  layer  cor¬ 
rection  will  satisfy  the  linear  system 


Ho"  *  'Vo  ' 


ht  -v o 


where  (F^,G^)  ”  ( fj (1 , 0) ,g2 ( 1, 0) )  .  The  decaying  solution  is 

-v 

Vq (c )  -  e  PjVgfO)  , 

-1  ( 50) 

nQ(o)  -  FjG^  vQ(o)  > 

where  matrix  P^  projects  onto  the  constant  n  k  dimensional 
unstable  eiqenspace  T  of  G^.  The  limiting  terminal  conditions 
take  the  form 

B(Xq(1)  ,G(X0(1)  ,1)  +  l?1v0(0),0)  -  0  .  (51) 

Assuming  that  we  can  solve  (perhaps  nonuniquely)  n  -  k  of 
these  equations  for 

P1v0(0)  *  6 (xQ Cl) >  ,  (52) 

there  will  remain  r  -  n  +  k  terminal  conditions 

T (XQ(1) )  5  C B ( Xq ( 1 ) , G ( Xg ( 1 ) ,1)  ♦  «(X0(1),0)}“  -  0  (53) 

for  Xq,  as  in  (25).  Thus,  any  limit  of  a  solution  of  the  form 
(44)  must  be  determined  by 

XQ  -  F(X0,t>  , 

(54) 

♦  <X0 CO) )  -  0  .  T(X0(1))  -  0  . 


Assuming  this  reduced  problem  is  solvable,  we  should  be  able 
to  construct  higher  order  approximations  to  solutions  without 
unusual  difficulty. 

we  wish  to  point  out  that  the  initial  impulse  i  uq(t)  « 

1  Got/c 

-  e  PqUq(O)  in  the  representation  (44)  of  the  fast  vector 
y  behaves  like  a  multiple  of  a  matrix  delta  function  in  the 
limit  c  ♦  0.  It  forces  the  rapid  transfer  (i.e.,  the  non- 
uniform  convergence)  of  the  slow  vector  x  from  x(0,0)  • 

X0!0.  +  =  XQ(0)  ♦  F2G0lpO;iO<0>  t0  V0>  in  the  * 

sional  range  1  or  the  stable  eigenspace  of  Gg.  This  is  in 
contrast  to  the  initial  Heaviside  jump  in  the  y  vector  and  the 
uniform  convergence  of  the  x  vector  for  asymptotic  solutions 
of  the  form  (4). 
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